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Abstract 

The physical sense and the properties of a group screening parameters which are determined 
in the previous papers (Part 1 and Part 2) for single- and two-component systems is discussed in 
this paper. On the base of data from the mentioned papers are determined two new characteristic 
lengths which complete a new hierarchy system of screening lengths in plasma. It was shown 
that the methods developed in Part 1 and Part 2 generates results which are applicable to the 
strongly non-ideal gaseous and dusty plasmas, and manifest a very good agreement with the existing 
experimental data. 
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I. INTRODUCTION 

In the previous papers , here Part 1 and Part 2, are developed the methods for 
describing of the electrostatic screening in single- and two-component systems (gas of elec- 
trons on a positive charged background, some electron-ion and dusty plasmas etc.). These 
methods generate a group of the new screening parameters which characterize the consid- 
ered systems. In this paper the physical sense of these parameters is discussed, and the 
reinterpretation of some known screening parameters (Landau radius, non-ideality parame- 
ters etc.) is given. On the base of results from Part 1 and Part 2 in this paper are obtained 
two new screening parameters which provide the possibility to form a new hierarchy system 
of the characteristic screening lengths and to compare the obtained results with the existing 
experimental data. 

The material presented in this paper is distributed in four Sections and one Appendix. 
In Section [TT] are given the some quantities and relations from Part 1 and Part 2 which 
should make easier reading of this paper; in Sections IIIII and IIVI are considered " small" 
characteristic lengths and connected with them non-ideality parameters for single- and two- 
component systems. Besides, in the same Sections are introduced "medium" and "large" 
characteristic lengths for the considered systems. The obtained results are compared with 
the existing experimental data in Section IIVI The conclusions of this paper are given in 
Section |Vl Finally, one important example of the application of obtained results, related to 
the systems with more than two components, is considered in Appendix |Al 



II. REMARKS: THE MAIN QUANTITIES AND RELATIONS 
A. Single-component systems 

The density, temperature and the charge of free particles in the initial (homogeneous) 
system: N, T and Ze, where Z = ±1, ±2..., and e - the modulus of the electron charge; 
The screening constant: 

[dfi/dN\ ' ^ ' 

where dfi/ dN = dfi{N, T) / dN, fi{N, T) is the chemical potential, and r^j- the corresponding 



characteristic length; 
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The charge density in the accessory (non-homogeneous) system with the probe particle: 



p{r)^ Pb + Ze-n{r), (2.2) 

where pb = —Ze ■ N is the background charge density, and n{r) - the free particle density 
at the distance r from the origin of the chosen reference frame; 
In the region r > Tq we have that 

p(r) = -ZeN ■ ro • exp(/«ro) • (2.3) 

r 

The radius ro and the parameters js{x) and 7k (x): 

ro = • 7s(a;), tq = r^'j^{x), (2.4) 



7,(a;) = [(1 + x')'3 - l]/x, j,{x) = (1 + x'^)^ - 1, (2.5) 



where is Wigncr-Seitz's radius (the radius of the sphere with the volume 1/N); 
The condition of electro-neutrality of the accessory system: 



oo 

Ze + 





J p{r) ■ inr'^dr = 0, (2.7) 



where Ze is the charge of the fixed probe particle. 



B. Two-component systems 

The density, temperature and the charge of ions in the initial (homogeneous) system: N^, 
Ti and Z^e, where Z^ — 1, 2...; 

The density, temperature and the charge of electrons in the same system: N^, Tg and — e 
or Zee, where Ze — —1; 

The ion and electron screening constants and the characteristic lengths: 

_ 1 r ^- 

K'i,e = = «0;i,e " (1 - «) 2 , Ko;i,e = 



^T^jZj^eey 
dpie/dNi 



(2.8) 
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where d/ii^e/dNi^e = diii^e{Ni,e,Ti^e)/dNi^e, IJ'i{Ni,Ti) and iie{Ne,Te) are the ion and electron 
chemical potentials, a is the electron-ion correlation coefficient given by expression 

a — 1 ^—^ (2.9) 

(1 + Xs) exp {-Xs) - (1 - Xs) exp (xs) ' 

and Xg — Ko-ef s;i'-i 

The charge densities in the accessory (non-homogeneous) systems with the probe parti- 
cles: 

p^'''\r) = Zie ■ nf''\r) - e ■ nf'''\r), (2.10) 

where nf'^\r) and ne'^\r) are the ion and electron densities at the distance r from the 
origin of the chosen reference frame; upper indexes {i, e) pointed the considered case: (i) - 
the probe particle charge is equal to ZiC, (e) - the probe particle charge is equal to — e; 
In the region r > rs;i,e we have that 

p(^'^)(r) = Ze,e(l - a) ■ iV,,ero,,e • exp(«:,ero;,e) • ^^P(~^'--^\ (2.11) 

where — —1. 

The radii ro;i,e in the cases (i) and (e) and the parameters 7s(a;i,e) ^-nd jK{xi,e)' 

Xi,e = l^i^eTs-M = rs;i,e = (^^) (2-13) 

where rs;i,e are Wigner-Seitz's radii (the radii of the spheres with the volumes 1/Ni and 

1/A^e); 

The electro-neutrality conditions of electro-neutrality of the accessory systems: 

oo 

Zi,ee + j p^''^\r) ■ Airr^dr = 0, (2.14) 



where Zj gC is the charge of the corresponding probe particle. 

III. THE SCREENING PARAMETERS OF THE SINGLE-COMPONENT SYS- 
TEMS 



A. "Small" characteristic length ro and the non-ideality parameters 7s ^ 

The connection of tq with Landau radius r^. As it was announced in Part 1, 
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all characteristic lengths of the developed method appear already in the case of single- 
component systems. The first of them is "small" characteristic length tq, given by Eqs. fl2.4p - 
(12. 6p . which is interpreted as a radius of the sphere centered on the probe particle and 
classically forbidden for other free particles. Because of that it was useful to compare the 
radius tq with a known characteristic length which has similar physical sense. Here, we keep 
in mind Landau's radius 



which is used in the case of the classical systems (see e.g. [3j). In accordance with Eqs. (12. 4p - 
(12. 6p we have that 

ro = ^K^rl ■ (1 + 0(x3)) , X < 1, (3.2) 
where in the classical case {dfj,/dN = kT/N) the relation 

-K^rl = ^—^ 3.3 
is valid. So, in this case the radii Tq and are connected by relation 

lim — = 1. (3.4) 

x^orL 

It means that represents an approximation od the characteristic length tq which is ap- 
plicable in the region of small x. Consequently, the parameter tq can be treated as the 
generalization of Landau's characteristic length which introduces into consideration from 
physical reasons. Let us emphasize that, contrary to ri which is principally unlimited, the 
radius ro < for any Kr^ > 0, in accordance with the conditions in real physical systems 
require, and that Wigner-Seitz's radius = lim^.^oo'"o- 

The connection of 7^ ^ with non-ideality parameters F and 7. From Eqs. (12. 4p - 
(ES]), ([32]) and follow the expressions 

Is = — , 1k = — , (3.5) 

for the coefficients 7^ and 7^. In connection with these expressions in the classical case 
one should keep in mind the relation (13.40 and the fact that = r^, where r^y is the 
corresponding Debye's radius. With respect to this it is useful to compare Eqs. (13.50 with 
the expressions for known quantities 

Ts kT-Ts'' Tf) kT ■ rn^ 



which usually use as the non-ideality parameters for the classical systems. Usually the 
parameters F and 7, similarly to Landau's radius ri, introduce from some physical reasons 
(see for example sl, Q]). From Eqs. fl2.4l) - fl2.6l) . (13.51) and (13.61) it follows that in the classical 

case 

lim 5 = 1, lim ^ = 1, (3.7) 

x^Q F x^o 7 

where x = nr^. It means that the parameters F and 7 appear here as approximative 
values of the coefficients 7^ and 7^ in the classical case in the region x ^ 1. Consequently, 
the parameters 7^ and 7^ can he treated as the generalized non-ideality parameters of single- 
component systems with any N, T and Z which make possible their non-relativistic treatment. 

From Eqs. (12.51) and (12. 6p it follows that the coefficients 7^ and 7^ represent the functions 
of the parameter r^/r^^ which is closely connected with another known classical quantity. 
Namely, in the classical case (r^ = td) we have that {rg/ri))''^ = n^,, where ud is so called 
Debye's number, i.e. the mean number of the free particles in the sphere with the radius r^. 

The behavior of the coefficients 7^ and 7^ and the parameters F and 7 as functions of 
is illustrated in Fig. [TJ This figure shows that in the region Krg < 0.5 the values of 
parameters F and 7 are practically same as the values of 7^ and 7^ and. They values left 
very closed up to the value = 1. 

B. "Medium" characteristic length Vc and the radius 

The charges Qin,out{f)- The electro-neutrality condition Eq. (12.141) can be presented in 
the form: Qin{r) + Qout{r) = 0, where < r < 00, and the quantities Qin,out{r) are given by 
expressions 

r 00 

QUr) = ^^ + j Pir'H^r'^dr', Qoutir) = J p{r')ATXT'^dr' , (3.8) 

r 

with the charge density p(r) defined by Eq. (12. 2p . One can see that Qmir) and Qouti.^) 
represent the total charge inside the sphere with radius r, centered in the point O, and the 
total charge of the rest of the space, respectively. 

Probably, the quantities Qin,out{r) could have different applications. So, the fact that in 
the considered case the total average charge of the probe particle self-sphere, i.e. Qin{r = r^), 
is strictly equal to the average charge of the free particles inside this sphere, can be of interest. 




FIG. 1: The behavior of the parameters 7s = ro/r^ and T = r^/r^, and 7^ = r^jr^ and 7 = rL/r£), 
defined by Eqs. p.4|) - (l2.6p and ()3.6p . as functions of x, in the classical case. 

However, in the further consideration the quantity Qin{r) will play only an accessory role. 
Namely, it will be used only the fact that Qin{r), obtained by means Eqs. (12. 2p and (13.81) . 
has the physical sense in the whole space (including the region r < r^), and can characterize 
the degree of the probe particle charge screening. 

The connection of Vc with the radius r^. In the case of single-component systems 
the screening constant k, defined by Eq. (12. ip . represents one of the main parameters in 
both methods - the method developed in this work and Debye-Hiickel's (DH) method. Con- 
sequently, the characteristic lengths = 1/k has to be treated in the similar way. 

Within DH method it is usual to interpret r^, which is equal to Debye's radius vd in the 
classical the screening radius. That is based on the fact that is a distance from 

the probe particle where the DH electrostatic potential (described in Part 1) becomes less 
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than Coulomb potential of this particle by the factor e~^. 

Here we will consider the role of radius from other aspect in order to clarify its real 
physical sense in both methods and determine the region of its applicability. For this purpose, 
we will introduce the radial charge density P(r) = 47rr^-p(r), where p(r) is given by Eq. (12 ■21) . 
and the characteristic length defined by relation 

|P(r,)| = max |P(r)|. (3.9) 

0<r<oo 

The behavior of P(r) for several values of x = nvg is shown in Fig. [21 This figure illustrates 
the fact that for /tr^ < 7^ the point r = there is in the region r > Tq, while for /tr^ > 7^ 
we have that = tq. Since in the region r > tq the charge density p(r) ~ exp{—Kr)/r, the 
parameter for Kr^ < 75 represents the root of equation 

dP{r) 



dr 



0, (3.10) 



where P(r) ~ r ■ exp(— /tr). Consequently, in this region the parameter is equal to r^. It 
means that within the method developed in this work we have the relations 

{r^, < KTs < 7^, 
, (3.11) 
ro, 73 < KTs < oo, 

which determine the real physical sense of the radius r^. Namely, from Eq. (13.111) it follows 
that has the physical sense only in the region /tr^ < 7^ where = Tc, while in the 
region nVg > 7^ the length loses any connection with the charge distribution in the probe 
particle neighborhood and, consequently, loses any physical sense. Let us emphasize that on 
the base of above mentioned one can conclude that it is always > tq. 

Then, we will consider the problem of screening of the probe particle as a problem of 
compensation of its charge Ze within the sphere of a radius r centered at the probe particle. 
For that purpose we will draw attention to the fact that uncompensated part of this charge 
in the case of such a sphere is equal to the charge Qin{r). Keeping this in mind, we compared 
the charge QmiTu) with the charge Ze in the region r^t > rg. Using the relations (13.81) and 
(12. 2p it could be shown that in this whole region QiniTu) > 0.735 • Ze. From here it follows 
that the radius r^, as well as the radius rc, cannot be treated as a characteristic length of 
full screening (neutrality) of the probe particle charge. Because of that, such a characteristic 
length is determined here in another way. 



8 




FIG. 2: The behavior of the radial charge density P(r) = 4TTr^p{r) for the cases: Kr^ = 1, when 
kvq < 1; Krs = 75, when kvq = 1; kVs = 3, when kvq > 1. The points where P{r) has the maximal 
value are denoted by rc- 

C. "Large" characteristic length r„ 

The charge Qin{r) and the quantity z/(r). In accordance with Eq. (I2.2p and the 
condition Eq. (12.71) . the charge Qin{r) can be presented in the form 

Qin{r) = -Qoutir) = Ze ■ u{r), (3.12) 
where the quantity z/(r) is given by the expression 

oo 

z/(r) = j [N - n{r')] ■ Anr'^dr', (3.13) 

r 

and n{r) is the free particle density in the considered system with the probe particle. 



Although the expressions (13.131) for z/(r) can be determined for any r > 0, this quantity 
has a special physical meaning in the region r > r^. Namely, it can be shown that in this 
region i/(r) represents the mean number of particles whose surplus into the mentioned sphere 
of radius r and corresponding deficit in the rest of the space causes the deviations of Qin{r) 
and Qoutij") from zero. Because of that only the quantities Qmij) and z/(r) in the region 
Ts < r < oo will be needed in this paper. From this reason in further considerations we will 
use the expression for z/(r) which is applicable in the region tq < r < oo, since it is always 
Tq < rg. In accordance with Eq. (20) for n{r) from Part 1 we have that 

z/(r) = (1 + Kr) exp (— Kr) ■ xi^)^ ^ > "'"o? (3-14) 

where x = kTs, and x{^) is defined by Eq (30) and illustrated by Fig. 3 from Part 1. 

Similarly to Qin{r), the quantity z/(r) could have also some different applications. For 
example, knowing of z/(r = r^) makes possible to estimate the density Np of pairs of particles 
with the inter-particle distance less then r^. Namely, in the binary approximation (not more 
than one particle inside the self-sphere of any particle in the initial system) it can be shown 
that: Np = (1/2) • • //(r^), where //(r^) is given by Eq. (I3.14p with r = r^. However, in 
further consideration the charge Qin{r) and the quantity z/(r) will play only an accessory 
role. 

The characteristic length r„ as the neutrality radius. In order to determine the 
required "large" characteristic length one should consider the charge Qin{f) = Ze ■ z/(r) in 
the region r > r^, where z/(r) is given by (I3.14p and has the sense which is described above. 
Let us remind that in this region the charge Ze of the probe particle is already completely 
compensated by background charge of the probe particle self-sphere. Consequently, in the 
region r > the charge Qj„(r) represents a quantitative characteristic of deviation of the 
neutrality of the sphere with radius r, centered at the prove particle, which is exceptionally 
caused by the charge of z/(r) particles which are enter in this sphere from the rest of space. 

From (I3.12P and (13.140 it could be seen that |(5m('")| in the region r > almost expo- 
nentially decreases from the maximum value \Qin{fs)\ down to zero, with the increasing of 
r. Consequently, as the requested screening length, denoted here as r„, can be taken the 
root of equation 




= e 



-1 



(3.15) 
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where Qmij) is given by fl3.12p and (13.141) . In accordance with this r„ can be treated as the 
neutrahty radius. Here, r„ is the radius of such a sphere centered at the probe particle for 
which charge exchange with the rest of the space becomes practically negligible in comparison 
with the similar exchange in the case of probe particle self-sphere, which is illustrated by 
Fig. [31 In the case of the initial system (see Part 1) the radius r„ can be interpreted as a 
radius of minimal sphere which can be considered as practically neutral one. 




FIG. 3: The visualization of the quantity v{r) = Qin{r)/Ze for r = rs and r = r, 



From (I3.15P and (13.141) it follows that the characteristic length r„ represents the root of 
the equation: 

1 + K,r 



Kr, 



exp[— /t(r — r^)] = e ^, (3.16) 
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which can be determined only numerically. Here, it is presented in two equivalent form, 
namely 

rn = rs-'ns, r„. = r^-T]^, (3.17) 
where the coefficients ris and 77^ are taken in the form 

r]s = [l + x + g{x)]/x, 1]^ = 1 + X + g{x), (3.18) 
since the member g{x) can be very well approximated by means of the easy expressions 

f 1.14619 -x + 0.43688-x2 < x < 1, 
g{x) = { ~ ~ (3.19) 

I 1/x - O.S/x^ + 0.08307/x^ 1 < x < 00, 

where at the point x = 1 both expressions give the same value g{l) = 0.58307. For the 
coefficients 77^ and r^^^ we obtain then following approximate expressions 
2.14619 



+ 0.43688 ■ X, < X < 1, 

X 



= \ ^ 1 / 1 0.5 0.08307\ ^ (3-20) 

1 + - 1 + - + — , 1 < X < 00, 

\ X x"' x-^ ) 



X \ X X 

2.14619 + 0.43688 -x^, < x < 1, 

riK= { ^ 1 / 0.5 0.08307\ ^ (3.21) 

X + 1 + -1 \ — , 1 < X < 00, 

X \ X x^ ) 

where x = kt^. The behavior of r^s and r/^ obtained numerically from basic equation fl3.16p 
and by approximative expressions (13.181) and (13.191) . is illustrated in Fig. HI This figure 
shows that the mentioned approximative expressions give very good results. From (I3.20p 
it follows that: limr/^ = 00 and lim r/g = 1. Consequently, on the base of (I3.17P we have 

a;— >0 x— >oo 

it that: < r„ < 00 for any kTs > 0, and Wigner-Seitz's radius = limx^oofn ■ Then, 



from (I3.2ip it follows that: > 1 for < nvg < 00. On the base of ^3.11^ we have it that: 
> f'c for any > 0. 



IV. THE SCREENING PARAMETERS OF THE TWO-COMPONENT SYSTEMS 
A. "Small" characteristic lengths rojjg and non-ideality parameters js;i,e and 7K;i,e 



The connection of ro;j,e with Landau radii rj^.j e- In Part 2 it was shown that 
the simplifications which give the possibility to describe a two-component plasma in DH 

12 




FIG. 4: The behavior of the coefficients rjg = Vn/rg, rj^ = rn/r^. Full lines show ijs and t/^ 
numerically determined from equation (|3.16|) : dashed lines - the same quantities determined by 
approximate expressions (|3.20p and (|3.2ip . 

method by means of unique screening constant are unacceptable and that its ion and electron 
components have to be described by own screening constant and Kg given by Eqs. fl2.8p 
and fl2.9p . One can see that these screening constants are substantially different from the DH 
screening constant (see Part 2). In the general case 7^ Kg. Unique exception is the case 
of completely classical plasma with Zi = 1, and Tj = Tg, but even in this case the common 
value of ion and electron screening constants is significantly different from DH screening 
constant. 

The existence of two special screening constants Kj and that in the case of two- 

component plasma we have two groups of screening parameters, analogous to the screening 
parameters described in the previous Section, but for ion and electron components separately 
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(see Part 2). As first, we will consider "small" characteristic lengths ro-i and ro;e, which are 
analogous to the characteristic length tq in the single-component case, and have the sense 
of radii of the spheres centered on the probe particles which are classically forbidden for the 
ions in the case (i) and for electrons in the case (e). 

Keeping in mind Eqs. (12.51) and (12.121) . as well as the facts that (1 — a) = 1 in the case 
of weakly non-ideal plasma and that in the classical case dfii^e/dNi^e = kTi^e/Ni^e, we have 
the relation 

lim 1:^ = 1^ (4.1) 
where Xi^e = i^i,ef's;i,e, and TL-i and ri-e are known ion and electron Landau's radii, namely 

~ kTi ' ^^'^ ~ kTe ' 
One can see that rm and ri-e represent the approximation of ro-i and ro;e in the classical 
case in the regions Xi^e ^ 1- Consequently, in the case of two-component classical plasma 
the characteristic lengths ro;i and ro;e can be treated as the corresponding generalization of 
Landau's radii r^.j and r^.g. It is important that, contrary to ri-i^e which are principally 
unlimited, the radii VQ-i e < r^.j e for any Ki^e'^s;i,e > 0, and that Wigner-Seitz's radii r^.j g = 

linixi,e^oo 0;i,e- 

The connection of 'js;i,e and 'jK;i,e with non-ideality parameters Fj e and 7^ e- On 

the base of Eqs. (12.51) . (12.121) . (14.11) and (14. 2 p in the case of classical plasma, we can obtain 
the relations 

lim 7lji:l = i^ lini 2W ^ (4 3) 

Xi,e^O Lj^e Xi^s-^O 7j^g 

where Xi^e = i^i,ei^s;i,e, and the quantities Fj e and ji^e are often used the classical ion and 
electron non-ideality parameters, given by relations 

i-i,e- ,rp , li,e - , „ • \^-^) 

'^-'■2,6 ' s;i,e '^-'■?,e ' K;i,e 

Let us emphasize that the parameters Fjg and 'ji^e, similarly to F and 7 in the single- 
component also introduced from "some physical reasons" (see for example M- 
However, one can see that Fj e and 7j^e represent the approximations of the ion and electron 
non-ideality parameters ■ys;i,e and 7K;i,e in the region Xj^g ^ 1- Consequently, 'js;i,e and 
7K;i,e can be treated as the generalization of ion and electron non-ideality parameters of 
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two- component classical plasma. The behavior of 7s;j,e and 7K;i,e as functions of Ki^efs;i,e is 
similar to the behavior of analogous parameters in the case of single-component system (see 
Fig. [ID . 

Let us draw attention that, contrary to the single- component case, the coefficients 'js;i,e 
and 7K;i,e illustrate the physical inadequateness of Debye's parameters Td and no for two- 
component systems. Namely, in accordance with Eqs. fl2.12p and fl2.13p these coefficients 
represent the functions of the parameters rs-^i^e/^K;i,e- However, in two-component case we 
have that r^^.^^e 7^ Td and the quantity (r^;i,e/r«;.j_e)"^ 7^ ^d- 



B. "Medium" characteristic lengths rc;i,e 

The charges Q^inluti^)- Similarly to the single-component case, we can take the electro- 
neutrality condition fl2.14p in the form: Q'fn\r) + Qout{f) = 0, where < r < cxd, and the 
quantities Qllf\r) and Qout{r) are given by relations 

r 00 

Qt^^) = Z,,ee + / p(^'^)(^') • 47rr'2rfr', Q^fir) = [ p^'''\r') ■ A7Tr'''dr', (4.5) 







with the charge density p(*''^)(r) defined by Eq. fl2.10p . Similarly to the single-component 
case, Q'llf\r) is the total charge of the whole sphere with radius r, centered at the probe 
particle, and Qo^\r) is the total charge of the rest of space. 

In all further considerations it is needed to know the charges Q^ijf\r) only in the region 
'^s;i,e < r < infty. Keeping in mind Eq. (12. lip , as well as the fact that Qi^^'*(r) = —Q'tut{r) 
under the condition (12.140 . we have that 

Q^n\^) = ^i,e^ ■ (1 - «) ■ X{Xi^e) (1 + K.i,er) exp {-Ki^er) , T^-i^e < T < OO, (4.6) 

where Xi^e = ^i,efs;i,ei = —1 and x{^) is the same function as in Eq. (I3.14p . 

For the practical applications of described method it is important that the expression 
(12. 9p for the coefficient of the electron- ion correlation a, which figure not only in Eq. (14.60 . 
but in expressions for all relevant quantities, can be very well approximated by two simple 
expressions, namely 

« = 3^ ■ a^s, < < I (4.7) 
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where the right side represents the first member of the expansion of right side of Eq. fl2.9p 
in the series with the respect to Xg, and 

a ^ —xl ■ (l + —xi) , < X < 10^/2. (4.8) 
The behavior of the coefficient a, given by Eqs. (I2.9p . (14 .Zp and (14.80 is illustrated by Fig. O 




X =K„ r 

s 0;e s;i 

FIG. 5: The behavior of the parameter a, given by exact expression ()2.9p . and by approximative 
expressions (j4.7p and (j4.8p . 

The connection of Vc-i^e with the radii VK-i^e- In the two-component case the ion and 
electron screening constants Kj and He have the sense of basic parameters of the method 
developed in Part 2. Consequently, the corresponding radii r^. = l/ni and r^^ = l/^e also 
represent the screening parameters of this method. In order to clarify the real role of r^^ ^ 
we will introduce the characteristic lengths Vc-i and Vc-e which in the cases (i) and (e) have 
the similar sense as the radius Tc in the single-component case. 
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It can be shown that in the regions < Ki^efs;i,e < 1 the parameters Vc-i^e are the roots of 
equations 

' =0, (4.9) 



dr 

where p(*'^)(r) = 47rr^ ■ p*^*'^^(r), and p(*'^)(r) is given by Eq. (12. lip . Consequently, we have 
that in the region < K,i^ef's;i,e < 1 the relations 

are valid. Due to the behavior of p^^'^\r), which is described in Part 2, we have it that 
in the > 1 the parameter rc-i^e < ''"s;i,e- In all examined cases in the region 

1 < Hi^efs;i,e < oo we obtain that 

except of the points Xi^e = 7^^'^, where it is Vc-i^e = r^-i^e = ro.i,e- From just mentioned it 
follows that the radii r^. and r^^ have clear physical sense only in the region < Ki^efs;i,e < 1 
where it is defined by the relation (I4.10p . 

Then, similarly to the single-component case we have to examine the behavior of the 
charge Q^in\f = ^c;i,e) for Ki^e'rs;i,e < 1- Siucc in this case Q\lf\r) is given by Eq. (14.61) . we 
have that Qfjf\rc-,i^e) > 0.735(1 — a)x{xi^e) ■ Zi(,e. It means that Q'fn\'''^c;i,e)i in accordance 
with the behavior of x{^i,e)i is comparable with the probe particle charge Zj^eC, exclud- 
ing eventually the region Xj_e ^ 1- Consequently, excluding the cases when a is close to 
unity, the parameters r^-^i and rc-ecannot be treated as a characteristic length of full screening 
(neutrality) of the probe particle charges in the cases of (i) and (e). Because of that, such 
characteristics lengths are determined here in the similar way as in the single-component 
case. 

Finally, we wish to draw attention to the fact that in the weakly non-ideal plasma (1 — 
a = 1) the radii r^.j and r^je are very closed to "medium" characteristic lengths in the 
corresponding single-component systems (ion gas on the negative background and electron 
gas on the positive background). This fact has already been considered in connection with 
non-applicability of DH method in the case of two-component plasmas. Namely, in p, la, 



m 



where the case of the classical plasma with Zi = 1 and Tj = Tg were considered 



(see Part 1), the screening constants and radius close to Ki^e and rK;i,e = l/'ti.e were used, 
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instead of DH screening constant kd and radius td = 1/kd- The results obtained in this 
work justify such a choice. 



C. "Large" characteristic lengths r^ j and r„;e as the neutrality radii 

The expressions for r„.j and r„.e. Similarly to Part 1, we will introduce the screening 
lengths r„;j and r„;e which represent the roots of the equations 

^ ' - = e-\ (4.12) 



Qin iTs;i,e) 

where the charge Q'in\r) is given by Eq. (14. 6 p and one of Eqs. (12. 9p . (14.70 and (14.80 . From 
it follows that Eq. (14.121) can be taken in the form 

^^'^'■''^ . exp [-fi;j,e(r -r^.j^e)] = e"\ (4.13) 



1 ~1~ ^i^e^ s\i 

which is the same as the form of the corresponding equation for the radius r„ from previous 
Section. Consequently, we have that r„;j and r„;e are given by relations 



where the coefficients ris{x) and rj^^x) are given by Eqs. (I3.18p - (l3.2ip . On the base of these 
expressions we have it that the relations rn;i,e > rs-i,e and r^^j^e = hnixi.eiooo i^n;i,e, as well as 
''^n;i,e > fc;i,e^ are Valid in the whole regions < Ki,e'"s;i,e < cxd. 

The behavior of the coefficients ?7s(xi,e) and rj^^Xi^e) in a wide region of Xi^e = i<'i,ei^s;i,e, 
which are given by Eqs. (EJ]), (EHl), (127[3D . (1338|) -( 13:2T]1 . (ITOIl and (ITOil . is displayed in 
Figs. [6] and [71 These figures clearly illustrate the changes of the behavior of f]s{xi^e) and 
Vnixi^e) with the increasing of deviation of the screening constant /tj^g from their classical 
values in the Namely, Fig. [6] shows the behavior of ris{xi) and rj^lxi) in the 

region < = Ko,e''"s;e < 2 where the classical and general expressions give practically the 
same results, while Fig.[7]is related to the region 2 < < 15 where the classical and general 
expressions give very different results. 

The applications and comparison with the existing experimental data. As we 
have already mentioned, the necessity of interpretation of experimental data caused several 
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FIG. 6: The behavior of the coefficients rjg and r/^ in the region < = i^Q-efs;e ^ 2 where the 
classical and general expressions for these coefficients give practically the same results. The case 
of plasma with Zi = 1 and Tj = Te, when Ko;e^s;e = i^O;e^s;i = Xs is presented. The shift of the 
crossing point of the curves presented related to the point KQ-ei's;e = 1 is caused by the application 
of the general expressions for rjs and ry^ . 

attempts (see for example [lO, 11, 12|) to determine of the characteristic screening length 
Tscr for two-component plasma which was taken as: r^cr = ^c;D ' ^d-, where r^) is Debye's 
radius for two-component system (see Part 2). In order to compare the values of r^cr and 
r„ we take here Tscr in the form 

^scr ' ^ K,\ii (4.16) 

where the correction factor kc = kc-oi^n/i^K^i- We will have in mind that the classical plasmas 
with Zi = 1 and Tj = Tg has been considered in the above mentioned papers, and that in 
such plasmas: Vg-i = Vg-e, Hi = He, r^-i = r^^-e, Tn-i = Tn-e and -q^-i = -q^-e- The comparison 
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(cl) 



FIG. 7: The behavior of coefficients rjg and in the region Xg = '«0;e?'s;e > 2. The curves rj. 
and rjli^^ show the behavior of these coefficients obtained in the classical case. The real behavior 
of the same coefhcients obtained with the general expression (j2.14p for kq-^ is presented with: □- 



for Np = 10^^ cm 



-3 



for Ne = W"cm--\ A- for = W^^cm~-^, and y- for = W^cm 



-3 ; 



m 



the case rjg, and with: 0- for = W^^cm ^, •- for Ng = IQ^^cm ^, for iVg = lO^^cm ^, and *- 
for N(, = lO^^cm"^ in the case rj^- 



of our results with the existing experimental data is performed for several cases and is 
presented in Fig. [SI This figure shows the behavior of the parameter ?7«..,- determined by 

nnnnnn 

fl4.14p - fl4.15p and the correction coefficients kc taken from [10|, [111, Il3|, [IJ, Il5|, Il6| in the 
region 0.6 < Kr» < 1.0. One can see a good agreement with the correction coefficients 

nn 

obtained from several measurements of conductivity from [141. |16I|. 

In order to examine usage of r„.j e as a screening radius in the expressions of Spitzer's 
type for the conductivity of plasma, we performed a calculation of conductivity for a fully 
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n^(x^=0.5)=5.657 2^''(xrj'n^=1 

n^(x^=1.0)=0.707 

n^(x^=1.5)=0.210 




x=Kr 

I I s;i 



FIG. 8: The comparison of the coefficient r/^ with the correction coefficients kc = rscr/^K;i, where 
rscr is the effective screening length determined in several papers on the base of experimental data. 
The cases of completely classic al p lasmas with Zj = 1 and 1^= Tg are presented. The values of kc 



are shown with: -k and *- 



Q,Q,Vand o- 



With © are shown 



the values of k^ for the same A'^e and Te as in [14], but determined by means of expression for the 
plasma conductivity from [17|. The curves /cc; KN and kr;nT,F! show the behavior of kc determined 



by means of analytical expressions from [10| and 



11 



161], respectively. 



ionized plasma with Zi = 1 and Ti = Tc = T, where r„.j = r^.g. The cases of plasmas 
with Nc = 10^^cm~^ and lO^^cm"^ in the region lO'^K < T < 5 • lO'^K were considered. 
The conductivities were determined by means of corrected Spitzer's expression |l^ where 
Debye's radius ro for two-component plasma is replaced (in so called Coulomb logarithm) 
by screening radius r^cr- The calculation was performed for 

^scr ^Dy ^scr 

where kxN is the corrected factor from [l^, and Vgcr = '^n^i- The results are presented in 
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Figs. M and [TOl In the same figures the corresponding values of conductivity determined by 
improved RPA method which is applicable for dense non-ideal plasmas P, y, 1^ are also 
presented. 
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FIG. 9: The Spitzer's and RPA conductivity of plasma with Zi = 1, Tj 



T and 



10 (cm ). The values of RPA conductivity (•) are taken form [12]. Spitzer's conductivities are 



determined by means of expression from 



11 



IJ], with the corrected screening radius rgcr, and the 
corresponding calculations are performed for: r^cr = I'd - dotted curve; r^cr = ■ kxN - dashed 
curve; r^cr = Tn-i - full curve. Here is Debye's radius for two-component plasma (see Part 2), 



and kxN is the corrected factor from 



The Figs. [9] and [TO] show that Spitzer's conductivity with r^cr = fn;i apparently converge 
to RPA conductivity with the increasing of electron density A'^e from lO^^cm"^ to lO^^cm"^, 
while Spitzer's conductivity with r^cr = and r^cr = kxN ■ rn lie significantly above. One 
can see that use of r^cr = rn-i provides the complete agreement with the results of RPA 
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FIG. 10: Same as in Fig. H but for Ne = W^^{cm"^). 

calculations in the case A^^e = lO^^cm^^. It is important that this agreement becomes better 
when the non-ideality degree increases. 

V. CONCLUSIONS 

The new model method for describing of the electrostatic screening in single- and two- 
components systems (electron-ion plasmas, dusty plasmas, some electrolytes, etc.) developed 
in Part 1 and Part 2 of this work, generates a group of new screening parameters. Here, 
we keep in mind "small" characteristic length vq and the non-ideality parameters 7^ and 
7k in the single-component case, and the corresponding characteristic lengths ro;i,e and the 
non-ideality parameters ■ys;i,e and 7K;i,e in the two-component case. 

In connection with the mentioned screening parameters is established that tq and ro;j,e 
represent the generalization of classical Landau's radii tl and rL-i^e, and ■Js,k, ls;i,e and 
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7K;j,e - the generalization of known classical non-ideality parameters F and 7 in the single- 
component case, and Fj^e and 7i,e in the two-component case. 

Apart of that, in this paper are introduced into consideration "medium" and "large" 
characteristic lengths and r„ in the single-component case, and rc;j,e and r„.j g in the two- 
component case. The behavior of these radii is examined in the whole regions < < 00 
and < Ki^efs;i,e < oo, whcrc K and Ki^e are the corresponding screening constants, and and 
rs;i,e- the corresponding Wigner-Seitz's radii. It was found that the considered characteristic 
lengths satisfy the relations 

ro < < r„, = lim tq = lim r„, (5.1) 

x—*oo x— »oo 



To-, 



lim 



lim r„ 



(5.2) 



ro < rc < r„, ro;i,e < rc;i,e < rn;i,e, (5.3) 

where x — kTs and xi^e — >^i,efs;i,e- These relations establish the two hierarchy systems of 
the characteristic lengths, and causes a redefinition of Wigner-Seitz's radii as a boundary 
screening lengths. 

Then, it was found that the radius = I/k in the single-component case has the sense 
only in the region of x < 7^/^, where r^, — Tc, and the radii rK;i,e — l/'^i.e have the sense only 
in the regions Xi^g < 1, where r^-i^e — fc;i,e in the two-component case. 

The results of application of the characteristic length r„;i as the neutrality radius were 
compared in this paper with existing experimental data in the cases of the classical plasmas 
with Zi = 1 and Tj = T^. It was found their very good agreement. 

Finally, we wish to draw attention that developed method is suitable for some astrophys- 
ical applications. Here we keep in mind that in outer shells of stars the physical conditions 
change from those which correspond to the rare, practically ideal plasma, to those which 
correspond to extremely dense non-ideal one. However, the method presented gives a pos- 
sibility to describe the electrostatic screening of all such outer shells in the same way, by 
means of the obtained screening characteristics. 
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APPENDIX A: THE LOWERING OF THE ATOMIC IONIZATION POTEN- 
TIAL IN PLASMA 

The method developed in Part 1 and Part 2 serve for describing of the electrostatic 
screening in the single- and two-component systems which contain only charged particles. 
However, as it is well known, the presence of the neutral component in plasma can often 
be neglected from the aspect of inner-plasma screening. It gives the possibility to apply 
the results obtained in this work on the problem of the lowering of the atomic ionization 
potential in plasmas with the neutral component (see Part 1). In such a way we will be able 
to compare the results obtained within the developed and DH method. The results of this 
comparison will give another important example of the inapplicability of DH method. 

In this context we will remind that in Part 2 the potential energies f/(*) and f/(^) of 
the ion and electron in plasma were determined. In the case = 1 and Tj = Te, when 
ion potential energy was compared with the corresponding DH ion 
potential energy U^jj' . Then, we will denote with /\Iat and /\Iat-D the atomic ionization 
potential determined within the developed and DH method, and take into account that 
A/at ~ and Alat-o ~ where the proportionality coefficients are equal or at least 
very closed. From here it follows the relation 

A/at 



(Al) 



The behavior of the right side of this relation is shown in Fig. 5 of Part 2. Accordingly to 



this figure we can conclude that Alat/ Alat-n < 0.8 in the who 
important to noted this fact, since the handbooks (see also 3 



e region < HiVs-i < oo. It is 
), often used in laboratories, 



recommend just DH lowering of the atom ionization potential Alat-^o, which would not be 
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taken as one of screening characteristics of electron-ion plasmas. 
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